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ABSTRACT 

o : 

S-j ' We present a search for the near infrared spectroscopic signature of the close orbit- 

ryj , ing extrasolar giant planet HD 75289b. We obtained ~230 spectra in the wavelength range 

2.18 \im- 2.19 \xm using the Phoenix spectrograph at Gemini South. By considering the di- 
rect spectrum, derived from irradiated model atmospheres, we search for the absorption profile 

■ signature present in the combined star and planet light. Since the planetary spectrum is sepa- 

■ rated from the stellar spectrum at most phases, we apply a phase dependent orbital model and 
CSl | tomographic techniques to search for absorption signatures. 

r**"- . Because the absorption signature lies buried in the noise of a single exposure we apply 

04 ■ a multiline deconvolution to the spectral lines available in order to boost the effective S/N 

ratio of the data. The wavelength coverage of 80 A is expected to contain ~ 100 planetary 
l/" - ) ■ lines, enabling a mean line with S/N ratio of 800 to be achieved after deconvolution. We are 

f^) ' nevertheless unable to detect the presence of the planet in the data and carry out further sim- 

C*"** . ulations to show that broader wavelength coverage should enable a planet like HD 75289b to 

be detected with 99.9 per cent (4 a) confidence. We investigate the sensitivity of our method 
and estimate detection tolerances for mismatches between observed and model planetary at- 
•i— I . mospheres. 

■ Key words: Line: profiles - Methods: data analysis - Techniques: spectroscopic - Stars: 
d ' late-type - Stars: individual: HD 75289 - Stars: planetary systems 



1 INTRODUCTION 

One of the great surprises in the search for worlds orbiting other 
stars was the discovery of a class of close orbiting Extrasolar Gi- 
ant Planets (CEGPs). The existence of such objects did not fit with 
the formation scenarios used to explain the arrangement of planets 
in our own solar system. Even before the first extrasolar gaseous 
giants were disco vered, core ac cretion models predicted that they 
should only form (Pollack 1984) at distances similar to those of the 
giant outer planets dBosslll995l) . Here temperatures must be cool 
enough to enable the formation of an ice (water, methane and am- 
monia in solid or liquid form) core from the proto-planetary disk 
before gas accretion onto the core can then take place. 

With the announcement of the discovery of the first CEGP or- 
biting the F8 dwarf 51 Peg (Mayor & Qu elozHl995h at a distance 
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of only 0.05 AU came the realisation t hat existing theories of plan- 
etary formation needed to be re vised. i Guillot et al] d 19961) argued, 
contrary to conventional belief dBosslll995h . that gas giant plan- 
ets could survive in such close proximity, although their formation 
beyond the ice point was still required. Mechanisms to enable an 
exchange of angular momentum to take place in a proto-planetary 
disk are believed to be responsible for the spiralling into and even- 
tual halting of planets in a close orbit dLin et alj 19961) . 

Planetary atmospheric physics research in recent years has 
been motivated by the expectation that radiation processed by a 
large body close to its parent star should be detectable with modern 
instrumentation. The fi rst clear detection of sod ium in the atmo- 
sphere of HD 209458b dCharbonneau et ai]|2002h revealed a lower 
than expected concentration when compared with a cloudless plan- 
etary atmosphere. Several explanations for the lack of sodium were 
put forward, including the possible presence of a high cloud deck. 
No observations had been able to differentiate between the possi- 
bility of no clouds (combined with lower sodium abundance) until 
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recent Spitzer obse rvations of HP 209 458b with the low resolution 
spectrometer IRS dHouck et al]|2004h in dicated a spectral featur e 
consistent with a high silicate cloud deck dRichardson et al]|2007t) . 
The reflected light s pectro sc opic studies carried out 

by. 



Phase function models (i=90°) 



Ihe reflected light s pectro sc opic studies carri e d out 
i Collier Cameron et al J 1 19991; ICharbonneau et ail 1 19991 ; 



ICollier Cameron et alJ 20021 ; Leigh et al 



l2003iJlbh and 



- .... . - ..... _ _ _ - more 

recent results from MOST photometry ~ dRoweetal.ll2006l) place 
albedo upp e r limits of 0.1-0.25 on the atmospheres of CEGPs. 
lLeighetal.1 d2003bh has placed an upper limit of 0.12 on the 
geometric albedo of HD 75289b, suggesting that if clouds are 
present at all, they are highly non-reflective. These are somewhat 
lower than the solar system gas giants, Jupiter, Saturn, Uranus 
and Neptune, which p ossess geome tric albedos of 0.46, 0.39, 0.60 
and 0.58 respectively dKarkoschkall 1994b . These observed albedo 
limits for CEGP s rule out the high silicate cloud investigated in 
the models of ISudarskv et alJ (I20001). and c ontradict more recent 
Spitzer observations dRichardson et al.l 2007). 

A number of attempts to model the emergent spectrum of 
CEGPs have also been made in recent years. Su darskv et al.l d2003l) 
have calculated spectra for planets with a range of orbital sep- 
arations from the parent star, including the so called class of 
'hot roasters'. These authors also presented spectra in the 0.4- 
5 u.m region for specific systems. For the well studied transit- 
ing planet, HD 209458b, they found a phase averaged planet-star 
flux ratio of \o%io(F v /F t ) ~ —3.10. A re-eva l uation of the 
ISudarskv etai] d2003l) models by iRichardson et alJ d2003l) which 
include phase dependent effects have yielded 2.2 \xm flux ratios of 
up to login (-Fp ~ —2.89 for a cloudless atmosphere in the 
case of HD 209458b. Here, the incident radiation is completely ab- 
sorbed a nd re-emitted on the d ay-side of the planet. A more recent 
study bv lBarman et all d2005h (BHA05), which models day-night 
gradients, has yielded results consistent with previous studies. For 
no redistribution of heat (i.e. re-emission of radiation on the day 
side), they find logio(Fp/F») ~ —2.90, while uniform redistribu- 
tion of the incident radiation yields logio(F p /F*) ~ —3.05 at all 
phases for HD 209458b. 

Several pioneering attempts were made to detect CEGPs in the 
near i nfrared (NIR) from the ground in the combined star and plane t 



light dWiedemann et al.l200ll ; ILucas & Rochel2002l ; ISnellenl2005h . 
Success in the IR was achieved with the Spitzer Space Telescope, 
with detections of a reduction in t hermal emission duri ng secondary 
transits of HD 209458b (GOV) dbeming et alj|2005h and TrES-1 
(K0V) dCharbonneau et al.ll2005h . HD 209458b photometry indi- 
cates logio (Fp/F„) ~ —2.59 at 24 trm while a 2.3 \im ground 
based measurement yields a value of logio(F p /Fi.) ~ —3.00, al- 
beit w ith a large degree of uncertainty ( Snellen 2005 ; D eming et al.l 
l2006h . The TrES-1 photometry dCharbonneau et alJ l2005h indi- 

cated logi (Fp/F„) <■ 3.15 at 4.5 tim and logi (F p /F») ~ 

—2.66 at 8.3 urni. T hese observati ons are consis t ent w ith the 
ISudarskv etai] d2003h base model. iDeming et al] d2006l) found 
that the CEGP HD 189733b yielded an even higher contrast ra- 
tio (logio(Fp/F*) ~ —2.26 at 24 (j.m) due to the close orbit and 
later spectral type of the K2V star. Fig. 8 of iFortnev et al.l (2006) 
plots a comparison of observed and predicted star-planet contrast 
ratios covering the spectral region 3 u.m- 30 u.m. The HD 209458b 
and TrES-1 observations are consistent with the models, with a 
2.2 u_m signal strength of 0.1%, whereas scaling the HD 189733b 
signal suggests a signa l strength as high as 0.2%. Most recently, 
Harrington et al. ( 2006J) have used high S/N Spitzer photometric 
observations to directly measure the planetary flux of the CEGP v 
And at 24 u.m. They find a phase dependent flux which is consis- 
tent with the BHA05 models where heat is not significantly redis- 
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Figure 1. Phase function, g(ce, A), representing three different models with 
an inclination of i = 90°. Plotted are an empirical reflected light Venus 
lHiltonlll992l) phase function, a Lambert Sphere phase function and the 
2.2 ujti phase function derived from the atmospheric models of BHA05. 
In the phase range 0.28 ^ 4> 0-72 , the BHA05 model closely mimics 
the Venus phase function model. 



tributed throughout the atmosphere but re-radiated on the dayside 
of the planet HD 209458b. 

1.1 HD 75289 

A companion to t he main se q uence GO dwarf HD 75289 was 
first announced by lUdrv et a?] d2000h following precision radial 
velocity monitoring with the CORALIE spectrometer at the 1.2- 
m Eu ler Swiss telescope. Since its first discovery, IPepe et aU 
d2002l) have published a refinement of the system parameters 
in light of further observations. A comprehensive list of sys- 
tem par ameters have been estimated and tabulated bv lLeigh et al.l 
d2003bh . Subsequent to further refinement of parameter s, as a 
result of a longer timebase of observations, iButler et aD d2006l) 
give an ephemeris of = 2450829.872 + 3.509267F day. This 
result, with P = 3.509267 ± 0.000064 day, is marginally greate r 
than the P = 3.50 91 ± 0.0001 day estimate of lUdrv et al] (2000). 
IButler et all d2006I) es timate an orbital radius semi-major axis of 
a = 0.0482 AU while lUdrv etal] d2OO0l) estimate a = 0.0483 AU. 

We have carried out Monte-Carlo simulations to determine the 
most probable estimates of planetary mass, Mp, and orbital incli- 
nation, i. An estimate of the orbital inclination provides an es- 
timate for the most probable velocity amplitude of the planetary 
signal. As described in Sj4]we can use this value as a guide when 
searching for a planetary signal. The axial inclination of the sys- 
tem can simply be calculated from the equatorial rotation veloc- 
ity, v eq , determined from estimates of the stellar radius and ro- 
tation period, while v sin i can be measured directly. We assume 
that the stellar equator and planetary orbit lie in the same plane in 
agreement with Rossiter-McLaugh l in effect measurements for the 



HD 2 09458 dBuiidv & Marcv|2000l ; IOueloz et alJ200d;IWinn et all 
l2005l ; IWittenmveret alJl2005l) and HP 189733 dWinn et al.ll2006t) 
systems. Monte-Carlo simulations using upd ated values d i scussed 
above and values tabulated in Table 1 of iLeigh et all (EfXBb) 
yield a most probable value of i — 67°. A planetary mass of 
M p = 0.51 ± 0.08Mj then follows from the known orbital am- 
plitude of the stellar reflex motion of K, — 54.9 ± 1.8 m s _1 . 
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Table 1. Journal of observations for 2006 January 15 & 16 Gemini South - Phoenix observations. Observations were made in an ABBA sequences. Each 
individual exposure was created from coadding shorter exposures. For HD 75063, the mean S/N ratio for each co-added group of spectra is given in column 
5. The combined HD 75063 spectra were used to remove the telluric features from each AB or BA combined pair of HD 75289 spectra post extraction (see 
Sj3 .41 . The mean S/N ratio after telluric lines and stellar lines were removed is given in column 5 for HD 75289. 



Object 


UT start 


Exp time per Number of ABBA 


Mean 


Comments 






exposure [s] 


sequences 


S/N 




2006 January 15 


HD 75289 


01:46:55 


120 


6 


401 


Seeing = 0.68" 


HD 75063 


02:50:32 


40 


2 


391 


G5V Telluric standard 


HD 75289 


03:18:48 


120 


12 


133 


Slit mis-alignment (counts decreasing) 


HD 75063 


05:22:47 


40 


3 


398 




HD 75289 


06:03:57 


120 


12 


380 


Low counts: slit re-alignment in middle of set. 


HD 75063 


08:46:42 


40 


3 


515 




2006 January 16 


HD 75063 


00:48:25 


40 


3 


414 




HD 75289 


01:14:02 


160 


10 


358 




HD 75063 


03:32:38 


40 


3 


240 




HD 75289 


03:59:03 


160 


11 


363 


Software crash at end of sequence: frames rejected 


HD 75063 


06:25:46 


40 


3 


385 




HD 75289 


03:48:51 


160 


11 


337 




HD 75063 


09:06:56 


40 


3 


405 





In light of observational and theoreti cal results, the adop ted 
theoretical 1.6-Rj radius determined by dLeigh et alj l2003ch is 
proba bly a significant over-estimate. A recent plot dBurrows et al.l 
2006) of radii determined from all observations of transiting 
CEGPs indicates 7? ~ 1 — lAR j. Based on the recen t theoretical 
evolution models from Fig. 3 of Burrows et alj ( feOOot) for CEGP 
radii, we estimate an approximate radius of R p ~ 1.2 Rj. 

In this work, we aim to detect the planetary signature of the 
close orbiting planet HD 75289b. In the NIR, rather than search- 
ing for an attenuated copy of the stellar spectrum, we can search 
for the direct spectrum emerging from the heated atmosphere 
wh ich is expected to po s sess temperatures in the range 1000-1200 
K dDeming et ID 120051 : ICharbonneau et al.l 120051 : iDeming et al.l 
2006). The analysis presented below relies on the use of the pre- 
dicted high density of absorption lines, due mainly to H2O, to 
search for the faint planetary signal. As such it has the potential 
to provide constraints on the reliability of model atmosphere opac- 
ities. In Sj2]we present the method used to detect a NIR planetary 
signal while Sj3] outlines the data reduction and analysis of data. 
We discuss our results in ij4]and carry out further simulations in Sj5] 
before discussing future prospects for this kind of survey in Sj6] 



In order to extract the planetary signal from a timeseries of 
spectra, we model the planetary signature as a phase dependent 
spectrum superimposed on an unvarying stellar spectrum. The time 
dependent variations of the planet orbiting the star are: (1) the in- 
stantaneous Doppler shift of the spectrum due to the relative posi- 
tion of the planet and (2) a phase dependent planet-star flux ratio, 
Fp/F*, which is dependent on atmospheric physics and heating 
due to the parent star. 

The ephemeris and velocity semi-amplitude, 
K p = VpSini — 137.6 km s" 1 (where i is the orbital incli- 
nation), enable the instantaneous velocity shift of the planetary 
spectrum relative to the stellar spectrum to be determined at any 
observation phase. Since the rotation velocity of HD 75289 is only 
3.8 km s _1 the planetary signature will be Doppler shifted clear of 
the stellar lines at all phases except those close to <j> = 0.0 & 0.5. 
A clean spectral separation will thus be present at most phases. 

Because we observe the combined light from star and planet, 
the light is dominated by the former, meaning that we must express 
the planetary signature as a fraction of the stellar signature. The 
planet-star flux ratio, e can be expressed as a function of orbital 
phase ((f>) and wavelength (A), such that 



2 DETECTING A NEAR INFRARED SPECTROSCOPIC 
SIGNATURE 

2.1 Phase Function 

The method of analysis uses a modification of the match ed fil- 
ter technique first presented in lCollier Cameron et alj d 19991) . This 
technique h as been refined and applied to several sy stems using 
optical data dCollier Came ron et al . 2002; [Leigh et al. 2003a. c) in- 
cluding HD 75289b dLeigh et alj |2003b) where an upper limit of 
login (-Fp/.F«) = —4.52 has been determined. Here we seek to ap- 
ply a similar method to NIR spectra. 



£ (a,A) 



/ P (t»,A) 
/.(A) 



en(A)ff(a,A). 



(1) 



The for m of this function is similar to that used by in optical 
studies (e.g. (Collier Camer on et al ] l2002l) ) but with the geomet- 
ric terms describing the albedo, planetary radius and orbital semi- 
major axis combined into a single function en(A). The value of 
eo(A) can be derived from model spectra at the desired value of 
A and represents the maximum planet-star flux ratio, observed at 
phase 0.5, when i — 90°. The phase angle a combines the orbital 
inclination, i, and orbital phase effects and is defined as 



cos(a) = —sin(i)cos((j>). 



(2) 
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The phase function g(a, A) may also either be empirically de- 
termined or modelled. 

For a cloud free model, lMarlevetal.ldl999l) showed that plan- 
ets reflect most efficiently shortward of ~ 6000 A, w here photons 
underg o Rayleigh scattering before being absorbed. [Seager et al.l 
(2000) found that the form of the phase function is strongly de- 
pendent on the particle size at 5500 A, with larger particles giving 
strong back scattering and albedos peaking at values similar to the 
Jovian planets in out own solar system. 

We have investigated the form of the phase function foi 
HD 209458 based on the results of BHA05. Fig. □ is a plot of 
g(a, A) for the case where i = 90° showing the close similar- 
ity of the Venus phase function and BHA05 phase functions for 
4> = 0.25 - 0.75. Since we are concerned with the detection of the 
planetary signal and not the characterisation (i.e. we are no attempt- 
ing to distinguish between phase function forms) in this paper, we 
have adopted the Venus phase function in subsequent analyses. This 
function has the advantage of being able to mimic inclination ef- 
fects through use of the a parameter, not possible with our em- 
pirical function, and is a valid approximation as our spectra were 
obtained with phases 0.275 ^ 4> ^ 0.72. 



3 OBSERVATIONS & DATA REDUCTION 

We present the obs ervations made with the Phoenix spectrograph 
( lHinkleetal1l2003h at Gemini South on 2006 January 15 & 16. 
Densely sampled timeseries were recorded using the 256 x 1024 
InSb Aladdin II array. In addition to the HD 75289 spectra, obser- 
vations of the bright Al III star, HD 75063, were made to enable 
careful monitoring of atmospheric absorption features. Because of 
vignetting and cosmetic issues, the useful region of spectrum was 
trimmed in the dispersion direction during processing to give a use- 
ful area of 256 x 850. The spectral range of 81.74 A covers the 
21807.32 A to 21889.06 A (2.18 urn -2.19 urn) region of the NIR 
at a spectral resolution of 56 800. With the 1024 pixel array, this 
gives a mean pixel resolution of 0.096 A per pixel which is equiva- 
lent to a 1.32 km s _1 velocity increment at the central wavelength. 
The journal of observations is recorded in Table Q] 



3.1 Detector 

The Aladdin array suffers from a number of significant cosmetic 
defects. These include a number of dust spots and elongated fea- 
tures. Hot pixels are also present and were flagged during the ex- 
traction process to ensure that they were not included. The left side 
of the detector appears to exhibit sensitivity which alternates be- 
tween adjacent rows. This stripe pattern produces a ripple at the 
~4% level but is not present on the right hand portion of the detec- 
tor. 

The characteristics of the Aladdin detector also necessitated 
an observing strategy to monitor the faint residual charge persis- 
tence signat ure common in this kind of detector when working in 
the infrared dSolomon| [l999). Observations were thus made by al- 
ternating the position of the star on the slit and thus the spectrum on 
the detector in an ABBA type pattern to allow any persistence sig- 
nature to be monitored. Each of the 4 sub-exposures in an ABBA 
sequence comprised a number of further sub-exposures. Short ex- 
posure times help to minimise dark current and may help to reduce 
persistence levels. For 40s exposures, 2x2 0s exposures were u sed. 
Similarly, 120s = 2x60s and 160s = 4x40s. lBeletic et al.ld2006f) has 
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Figure 2. Top: Aladdin II array master flatfield image showing bad and 
hot pixels, the large scale variations and the row to row pixel sensitiv- 
ity alternation on the left had part of the detector. The large scale varia- 
tions were removed post-extraction by fitting splines to the continuum. The 
greyscale runs from black=0.75 to white=1.25. Bottom: The mean profile of 
the largescale variations created after rotating the frame (to ensure the spec- 
tram ran parallel to the detector) and collapsing the spectrum. The spec- 
tram was then normalised by dividing by a straight line which was used to 
remove the variation in throughput of the spectrograph (from left to right). 
The plot shows that the mean variation of the large scale ripples is expected 
to be at the 5% level in the extracted spectra. 



however found that continual flushing of the array can increase the 
dark current floor, and increase the persistence current. 

We could not detect persistence of the HD 75289 spectra 
which typically peaked at 5000 e~s. The HD 75063 spectra how- 
ever yielded a maximum of 12500 e~s and did leave a faint per- 
sistence trace at a level of 1.6a of the sky background or 0.13% of 
the peak counts when switching between A and B positions. For 
HD 75063, telluric line s may shift wavelengt h throughout the night 
by a few tens of m s _1 dGrav & Browrw2006l) . but this effect trans- 
lates to a shift of l/130th of a pixel. Hence the persistence effect 
when coadding these spectra is negligible. The effect is more cru- 
cial, if significant variation in the position of features is expected, 
such as from an orbiting planet whose signature is Doppler shifted 
during its orbit. If a signal of similar relative strength to that in the 
HD 75063 spectra is seen in the HD 75289 spectra, the blurring ef- 
fect is still expected to be a second order phenomenon and we do 
not expect to detect such a signature. 
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3.2 Data Reduction 

Pixel to pixel variations were removed using flat-field exposures 
taken with an internal tungsten reference lamp. In order to create 
a reliable balance frame to remove the pixel sensitivity variations, 
we divided a Gaussian blurred (using a FWHM of 7 pixels) version 
of the master flat field image by the original master flat field image. 
Fig-IH(top) shows the Hatfield image normalised to unity. The large- 
scale variations at the 5 per cent level shown in Fig. [2] (bottom) are 
a consequence of broad sensitivity variations on the Aladdin II ar- 
ray. We found that these variations were dependent on the illumi- 
nation of the array, being both source and time dependent, result- 
ing in an inability to remove the pattern using standard flatfielding 
techniques. For example, we extracted the flat-field variation with a 
profile used for extraction of our object frames. The morphology of 
the flatfield ripples was found to differ from the continuum ripple 
in the extracted object frame. Similarly, since the ripple pattern was 
source dependent, we were unable to flatfield the spectra by divid- 
ing by the standard star, HD 75063. While the pattern for a given 
object remained stationary to first order throughout each night of 
observations, second order time dependent variation was also seen. 
Nodding the telescope (see below) between two positions on the slit 
resulted in a drift of the A and B positions of several pixels through- 
out the night. During extraction we therefore traced each spectrum 
independently. Since the ripple pattern described above appeared 
different for HD 75289b than for the flat field, we removed it by 
fitting splines to the continuum of the extracted spectra. 

The worst cosmic ray events were removed at the pre- 
extraction stage using the Starlink FIGARO routine BCLEAN 
dShortrideell 1993b , Instead of making use of the ABBA sequences 
to reject sky lines by extracting from the A-B and B-A differences 
(Joyce 1992), we found that the S/N ratios of the extracted spectra 
were optimised when the sky background was modelled by fitting 
polynomials of degree 3 to the pixels either side of the spectral 
profile at each X-position in each frame. An iterative fitting was 
used to reject deviant X-position sky fits, thereby rejecting any 
sky lines. Only a single sky line at ~ 21868 A was present in the 
data at a level of ~ 9 per cent in the worst spectra (on January 15) 
and 1 per cent in a typical spectrum.. The spectra were extracted 
using ECHO MOP, the echelle reduction package developed by 
I Mill si d 19941) . The spectra were extracted using ECHOMOP's 
im plemen t ation of the OPTIMAL extraction algorithm developed 
by iHornel {l986). ECHOMOP propagates error information based 
on photon statistics and readout noise throughout the extraction 
process. 



3.3 Wavelength calibration 

At the time of observations, no arc lamp giving sufficient lines to 
perform a calibration in the 2.18 |i.m-2.19 p.mwavelength range 
was available. We were thus unable to carry out a conventional cal- 
ibration using the positions of known emission lin es. Instead, we 
used a spectrum generated from a HITRAN line list (Rothm an et all 
2005) to identify corresponding features in an observed spectrum 
of the telluric standard star, HD 75063. We used the simple emis- 
sion line fitting routine EMLT in the Starlink FIGARO package to 
fit Gaussians to the inverted spectra in order to identify the posi- 
tions and width of the lines. The latter were used as uncertainties 
in the measurements. By fitting a cubic polynomial to the 12 iden- 
tified telluric lines (Fig. [3}, we obtained a wavelength calibration 
which was subsequently used in our analysis. The \ 2 of the fit us- 
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Figure 3. Residuals from wavelength calibration using the theoretical water 
vapour line positions measured from a HITRAN spectrum and the observed 
pixel positions of the corresponding features in the observed spectrum of 
our telluric standard star, HD 75063. The best fit was obtained using a cubic 
polynomial. 

ing a cubic polynomial provided improvements by a factor of 31 
and 1.8 when compared with a straight line and a quadratic fit re- 
spectively. The 0.023 A rms residual scatter in the fit corresponds 
to 0.32 km s" 1 at the centroidal wavelength of 21484 A, or 0.06 of 
a resolution element.. 

3.4 Extracting the planetary signal 

We first constructed a high S/N ratio master stellar spectrum tem- 
plate to accurately subtract the unshifted starlight from each ob- 
served spectrum. This has the additional benefit of removing (to 
first order) the telluric water vapour lines prevalent in this part 
of the spectrum. By monitoring the strengths of the telluric lines 
alone, we found, as expected, that they vary in strength throughout 
the night, generally being stronger when HD 75289 was observed 
at high airmass. The telluric lines did not all behave in this way 
however, with some lines weaker than the mean at times while oth- 
ers were stronger. Using a single master frame comprising both 
stellar and telluric lines left residuals of up to 5-a times the ex- 
pected noise level after subtraction of the starlight from each spec- 
trum. 

Instead we made individual use of the bright star, HD 75063, 
which contains only telluric lines in the observed spectral range 
of our data. This star was observed at regular intervals on each 
night of observations. A mean telluric template spectrum was cre- 
ated for each night of observations by aligning and coadding all 
the HD 75063 spectra observed throughout that night. The telluric 
spectrum was fitted to each HD 75289 spectrum by taking deriva- 
tives of the spectra and using splines to calculate the scale factor at 
points across the spectra. This process can account for lines which 
behave independently over the night (i.e. all telluric lines do not 
necessarily vary in strength by the same factor at any given time) 
and is described in ICollier Cameron et alj d2002l) (Appendix A). 
The procedure provides a spectrum which closely matches the tel- 
luric lines for each observed HD 75289 spectrum and is used to di- 
vide out these features, leaving only stellar and planetary features. 
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Figure 4. Model and observed spectra in the 2.18-2.19 |am span of the observations at a resolution of R ~ 56, 800. Top: Model planetary spectrum with tick 
marks below indicating the 8 strongest lines (see §3.5\ . Middle: Observed HD 75289 spectrum plus undetected planetary spectrum and telluric features (first 
observation from January 16). Tick marks below indicate the theoretical positions of Fe and Si stellar lines as given by VALD (see £13.51 . There are additional 
absorption features not listed in VALD at A = 21816.4 A, 21857.8 A and 21882.8 A. Bottom: HD 75063 Al III spectrum showing only telluric lines. The 
observed HD 75289 spectrum and model planet spectrum are offset by arbitrary units of 0.5 and 1.0 respectively for plotting purposes. 



A master stellar-plus-planet spectrum is then created for each 
telluric corrected HD 75289 observation, but excluding the current 
telluric corrected HD 75289 spectrum. Every stellar-plus-planet 
spectrum is then divided by an aligned and shifted version of the 
current master stellar-plus-planet spectrum in order to remove the 
stellar lines. Since a planetary signal is Doppler shifted according 
to phase relative to the stellar lines only the planetary signature 
should remain. The master stellar-plus-planet spectrum which is 
subtracted from each spectrum in turn contains a blurred copy of 
the planetary signal. This will thus attenuate the planetary signal in 
the frame on which subtraction is being performed, to some de- 
gree. The effect is worst at (f> = 0.25 & (j> = 0.75 where the 
planet moves very little in velocity space from spectrum to spec- 
trum. The attenuation correction applied to the matched travelling 
Gaussian used to model the planetary motion when searching for 
a signal is described in deta il in §6 and Appendices D2 & D3 of 
ICollier Cameron et alj f20o3) . 

We found however that further time varying systematic resid- 
uals remained in the spectra and believe that these remaining resid- 
uals are due to the time varying ripple which we were unable to 
remove from the spectra at the extraction stage. Since we remove 
this effect using splines, we expect time varying differences in the 
fits which do not exactly match the observed changes in the ripple 
pattern due to finite signal to noise in the data. We therefore imple- 
mented a method which removes the residual signal which we treat 
as pattern nois e. This method using princip al component analysis 
is described bv lCollier Cameron et alj ([2002) (Appendix B) and re- 
moves time varying patterns at fixed positions in the spectra. When 



the first few principal components are removed, this has little or 
no effect on any planetary signal since it changes position from 
spectrum to spectrum. We discuss the results using this algorithm 
further in section SB] 



3.5 Model spectrum 

The residual spectrum will contain a Doppler-shifted copy of the 
planetary spectrum, which at this stage is still buried in the noise. 
To reduce the effective noise, we model the planet's spectrum as 
the convolution of a Doppler-shifted mean line profile (which we 
wish to determine) and an absorption line pattern for an irradi- 
ated model atmosphere spectrum with parameters equal to those 
of HD 75289b. Our model spectra were generated using the cloud- 
free "rainout" method described in BHA05 which improves the ear- 
lier models of Barman et al. I feOOlh (based on the ' AMES-cond' 
models o f lAllardetal.[ ]2001)) by iteratively reducing the elemen- 
tal abundances involved in grain formation at a given layer and 
recomputing the chemical equilibrium at each new set of strat- 
ified elemental abundances. The resulting equilibrium chemistry 
and opacity sampling after rainout of species is fully self consistent, 
whereas earlier cond models simply excluded the grain opacities. 
Thus in the models used in this work, elements which are impor- 
tant in the atmospheres of cool stars such as Ti and V are signifi- 
cantly depleted leading to negligible concentrations of TiO and VO. 
Day-night temperature gradients were modelled under the assump- 
tion that concentric neighbouring rings of different temperature in- 
teract very little via radiative transfer processes. For a fuller de- 
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scription of the mod el opacities and setup see Allai'det al. (2001), 
Barm an et al, I< l200ll) and BHA05. The temperatures on the dayside 
of the planet lead to an atmospheric chemistry dominated by H2, 
He, H 2 and CO (BHA05). The dominant opacities in the 2.18- 
2.19 \im wavelength span of our data however are due to H2O. 
We derived the necessary list of line positions and depths for least 
squares deconvolution from a model spectrum before instrumental 
and stellar line broadening were added to the spectra and for the 
case when the planet is in conjunction with the star (i.e. at orbital 
phase 0.5). 

In Fig. [4] (top), we plot the model spectrum generated for 
HD 75289b. The spectrum has been convolved with a Gaussian 
to mimic the R ~ 56, 800 resolution of the observations and rota- 
tionally broadened to 1.482 km s~ under the assumption that the 
planet is tidally locked to the star (i.e. v sin i = 2nR p sin i/P). Fig. 
|4]also shows a typical HD 75289 spectrum (middle) and HD 75063 
spectrum (bottom). The tick marks below then HD 75289b spec- 
trum indicate the theoretical positions of 10 stellar lines (Fe and Si 
opacitie s) which are listed in the Vienna Atomic Line Database, 
VALD jKupkaetalJll999h . Additional opacities not included in 
VALD can also be seen (see Fig.|4]caption). The HD 75289 spectra 
are thus clearly dominated by telluric lines. 

The mean absorption profile is recovered through use of a least 
square s deconvolution method first demonstrated by Don ati et al.l 
d 19971) . A list of absorption lines is derived from the model at- 
mosphere (see above) and the depths, before any kind of broad- 
ening, are used to optimally weight each line in the observed 
spectrum. The deconvolution is effectively a sophisticated cross- 
correlation method which optimally coadds all the aligned absorp- 
tion profiles while removing any side lobes due to line blend- 
ing. The resulting profile is a mean absorption line with S/N ra- 
tio boosted, thereby improving the chance of detecting weak sig- 
nals which are dom inated by noise. The current version of the code 
jBarnes et al.ll998T) which propagates errors from the input spectra 
ha s been used extensively for reflected light searches in the opti cal 
bv lCollier Camer on et al. Ul999ll2002l) and lLeighetai] ( l2003al lb1). 
For HD 75289b, we expect ~98 planetary absorption features in 
the wa velength range of o ur data. The weighted mean normalised 
depth jBarnes et alj|l99§) of these lines relative to the planetary 
"continuum" is 0.096. The optimal nature of the deconvolution 
means that all 98 lines in the wavelength range of our observations 
contribute, but only 8 (indicated by tick marks below the model 
planet spectrum in Fig. |4) of the lines possess a normalised depth 
relative to the planetary continuum greater than the mean weighted 
value of 0.096. Of the 98 lines, 50 possess depths greater than 0.01 
of the normalised planetary continuum. Attenuation by a further 
factor in the combined star+planet spectrum is one of the main pa- 
rameters which we wish to determine. 

From photon statistics {\fW2 per spectrum, where N e = num- 
ber of electrons), we expect an mean S/N limit in our spectra of 
~ 269 for both nights of data. The mean S/N ratios, measured from 
flat regions in the normalised spectra, were 1991*55 ant ^ 268^77 
for the nights of January 15 & 16 respectively. The upper and lower 
limits represent the minimum and maximum recorded S/N ratio on 
each night. Slit alignment and tracking problems were the cause of 
low counts in a number of the January 15 spectra with a minimum 
S/N = 25 (readout noise = 40 e") . The gain in S/N of 3.4 for the 
deconvolved line profile yielded S/N ratios of 797 t.gg 2 over the 
two nights (i.e. minimum S/N ~1 12 and maximum S/N = 1 149). 



4 RESULTS 

4.1 Matched-filter analysis 

We model the time-dependent changes in Doppler shift and bright- 
ness, in the manner described in !}2] This matched-filter enables us 
to search for features in the timeseries spectra and is described in 
ICollier Cameron et alj J2OO2I) (appendix D). We asses the relative 
probabilities of the \ 2 fits to the data by varying e(A) and K p and 
plotting the improvement in \ 2 , which is normalised to the best- 
fitting model. To calibrate any candidate detection, we construct 
a simulated planet signal of known eo (A) that is added to the ex- 
tracted spectra prior to removal of stellar and telluric features and 
before least squares deconvolution is carried out. By ensuring the 
fake planet is recovered correctly by our procedures, we can be con- 
fident of a calibrated detection in the presence of a genuine planet 
signal. 

The significance of the result is assessed using bootstrap sta- 
tistical procedures based on the random re-ordering of the data in a 
way that scrambles phase s while preserving the effects of any cor- 
related systematic errors dCollier Cameron et al.l l2002). The order 
of the observations is randomised in a set of 3000 trials which will 
scramble any true planetary signal while enabling the data to retain 
the ability to produce spurious detections through the chance align- 
ment of systematic errors. The least squares estimate of logio e(A) 
and associated x 2 as a function of K p enable us to plot 68.4, 95.4, 
99.0 and 99.9 per cent bootstrap limits on the strength of the plan- 
etary signal. 



4.2 Searching for a planetary signal 

Fig. H] (left) shows the phased timeseries deconvolved spectra and 
the resulting x 2 map of logio e(A) vs K p . The effective wavelength 
after performing deconvolution is A = 21848 A. Fig. [5] (right ) 
shows the relative probability x 2 ma P °f the planet-star flux ra- 
tio logio e(A) vs K p . The darkest regions in the map repre- 
sent the greatest improvement in x 2 when fitting the model de- 
scribed in $4. 11 We can obtain an estimate of the upper limit 
(i.e. when i = 90°) of the velocity amplitude of the planet. With 
a = 0.0482 ± 0.0028 AU and P = 3.509267 ± 0.000064 d, we 
find v p = max(K p ) = 149.43 ± 8.68 km s _1 . With the orbital 
inclination of i — 67° reported in jyTjJ we estimate a most prob- 
able velocity amplitude of K p = 137.55 ± 7.99 km s" 1 . In Fig.f5] 
we plot the corresponding motion and recovered position for K p 
(dashed curves and vertical lines respectively) as a visual guide to 
the reader. 

It is possible to see the non-Gaussian systematics in the phase 
plot in Fig. [5] (top left). We believe that these features are the 
result of our inability to fully remove the varying telluric signa- 
tures and uncharacterised changes in the observed ripple function 
described in ^3.21 . To remove the most significant trends in the 
timeseries data at fixed positions in wavelength, we u sed princi- 
pal component analysis dCoilier Cameron et al.l (2002). Appendix 
B). The eigenvalues of the first two principal components indi- 
cated a strong contribution while a less strong contribution from 
the next five components was then followed by a linear decrease 
in significance. We thus investigated removal of the first two, and 
the first five components, with the resulting input timeseries for 
the latter case being significantly cleaner of systematic noise fea- 
tures. When we compare Fig. [5] (bottom right) with Fig. [5] (top 
right), we see how the K p — 157.8 km s _1 feature is reduced in 
significance, leaving several features of comparable strength and 
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Figure 5. Left: Phased timeseries of the deconvolved residual spectra of HD 75289b. The dashed sinusoidal curve represents the motion of a planetary signal 
based on the most probable velocity amplitude, K p = 137.55 km s —1 , estimated from empirically determined parameters ( i]4.2t and does not represent a 
detection with this amplitude. Right: Relative probability \ 2 ma P of planet-star flux ratio logio e(A) vs K p . Results are shown after removal of the first two 
(top) and first five (bottom) principal components from the timeseries spectra. The greyscale in the probability maps represents the probability relative to the 
best fitting model in the range white for to black for 1. Plotted are the 68.4, 95.4, 99 and 99.9 per cent (bottom to top) confidence levels for a planetary 
signal detection. The dashed vertical line represents the most probable velocity amplitude. The dark feature at K p = 157.8 km s _1 and log\a{e) = -2.2 is 
detected with < 95.4 per cent confidence. Note how removal of five principal components reduces the significance of this feature such that it becomes 
approximately equal in significance to the K p = 116.5 km s _1 candidate (bottom right). These features are close to the noise level and are probably due to 
incomplete removal of the fixed pattern seen in the timeseries spectra. 
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low statistical significance. The K p = 157.8 km s _1 feature only 
marginally lies in the range of possible values since the uncertainty 
on our v p = max(K p ) measurement yields an upper (1-cr) limit 
of 158.1 1 km s" 1 . A planet with this velocity amplitude would im- 
ply i ~ 86.5°, indicating a planet likely to show transits, an ob- 
servation not reported in the literature. The K v — 116.5 km s~ 
feature is of no greater significance and implies an orbital inclina- 
tion of i = 47.5° respectively. This is unlikely to give a detection 
given the relatively low maximum flux we wou ld receive under the 
reasonable assumption dHarrington eTaT1l2006h that heat is not ef- 
fectively redistributed and the planet is hottest on its dayside. All 
candidate features however can confidently be rejected on the ba- 
sis of the observed flux ratios for a number of systems (SjT) since 
the mean logio e(A) = -2.1 impl ies a planet significa ntly hotter than 
even HD 189733 (see Fig. 8 of lFortnev et"aH < l2006h ). 

Under the assumption that all candidate features are spurious, 
the overall level of noise in the data enables us to rule out the pres- 
ence of a planet at K v with logio £o(A) > —2.1 at the 2-<r level 
(Fig.0 top right). 



5 SIMULATING A PLANETARY SIGNAL 

In order to assess the feasibility of detecting a planetary signature in 
the NIR, we have carried out simulations using fake data based on 
the known physical parameters of HD 75289b. By considering the 
wavelength coverage possible with current NIR multiorder cross- 
dispe rsed echelle technology (i.e. IRCS/Subaru {Kobavashi et al.l 
l200d) and NIRSPEC/Keck 2 dMcLean et alj|l998h ). we estimate 
that approximately 50 spectra with R ~ 20000, a mean S/N ratio 
of 300 and encompassing a wavelength range of 1.90 \im to 2.45 
ixm c ould be obtained per 8 hrs of observations ( Kobavas hi et al.l 
2000). With a typical spectral coverage of 70% of this range our 
planetary model atmosphere for HD 75289 yields ~ 2300 useful 
lines with depths between 0.05 and 1.0 of the normalised contin- 
uum. 

We g enerated mod el spectra for a GOV star using ATLAS9 
models dKuruczl fl993h and a HITRAN model of the tel- 
luric lines dRothman et all l2005h . These spectra were com- 
bined with our model planetary spectrum assuming an or- 
bital inclination of i — 67°, a relative maximum planet-star 
flux ratio of logio eo = —3 and a velocity amplitude of 
K p — VpSini — 137.6 km s~\ Gaussian random noise was added 
to the spectra to simulate a S/N ratio of 300. We assumed observa- 
tions were taken over two nights with a ten night phase gap. This 
strategy optimises observations at phases close to <j> = 0.5, while 
minimising the timebase. 

We attempted to recover of the planetary signal in the same 
way as described in previous sections. The phased timeseries spec- 
tra are presented in Fig.|6](top left), and show a detection well clear 
of the 99 per cent confidence level (Fig. [6] top right). Since the 
semi-major axis of the orbit and the period are both known, we 



matched-filter analysis gives K p — 142.7 kms -1 leading to an 
inclination determination of 72° . Hence we overestimate the incli- 
nation by 5°. The detection limit at the simulated velocity of the 
planet is logio £o = —3.2. Equivalently, at the projected velocity 
of the planet, we expect to recover the signal at the limiting confi- 
dence for data with S/N ^ 200. 

We carried out a further test to estimate the feasibility of re- 
covering the planetary signal in the case where the line strengths 
are incorrect. We assumed that the positions of all lines in our input 



model planetary spectrum could be incorrect by some fraction. The 
depth of each line, but not the position, was modified by a Gaus- 
sian random fraction of the original strength. We then attempted 
to recover the planetary signal with the modified line list (used for 
deconvolution) and found that for data with S/N = 300, we still re- 
cover the planet with a 99.9 per cent confidence limit for strengths 
which have been modified by on average 50 per cent. This simula- 
tion essentially mimics the effects of non-optimal extraction of the 
lines using our deconvolution code. Our ability to recover the planet 
is more sensitive to incorrect positioning of lines however. We find 
that a clear 99.9 per cent confidence detection is achieved in the 
limiting case of randomly re-positioning 15 per cent of the lines be- 
fore deconvolution. The main reason these mismatches, which we 
discuss in the next section, is likely the uncertainties in opacities, 
opacity oscillator strengths and atmospheric physics of the planet. 
We note that in reality, removal of the telluric lines at the shorter 
wavelengths of the above range (i.e. 1.90 vim to 2.08 \im) does not 
affect our ability to make a clear detection. 

We also simulated recovery of the planet for a spectral reso- 
lution of R ~ 40000. The main advantage of high resolution data 
would be the increase in depth of the absorption lines relative to 
the continuum. The results are shown in Fig. [6] (bottom) and show 
the clear gain attained by doubling the resolution. Our limiting 99 
per cent confidence threshold in this instance would enable us to 
detect a planetary signal at a level of logio eo — —3.5. The planet 
is recovered with K p = 139.7 km s _1 , leading to an inclination 
determination of 69.2°. Since a narrower slit will be used to make 
observations at higher resolution, we may not observe this degree 
of improvement in reality. Our limiting 99.9 per cent confidence in- 
dicates we would still detect the planet with S/N = 95 or approxi- 
mately 1/4 the number of photons required in the R ~ 20000 limit- 
ing case. We discuss the prospect of obtaining spectra of this nature 
in section 36] 



6 DISCUSSION 

The methods employed in this paper provide a very encouraging 
indication that a planetary signal can potentially be extracted from 
high resolution timeseries spectra. The fact that we do not detect a 
signal can be attributed largely to systematic noise introduced from 
difficulties in consistently normalising each spectrum in our time- 
series. This arises largely because of the systematic and varying 
large scale flat fielding errors associated with the detector employed 
in the observations as discussed in 33.21 

We have shown that the feature detected in our spectra is un- 
likely to be a real planetary feature since its strength is much re- 
duced when additional principal components, which describe (sys- 
tematic) trends in the data, are removed. The eigenvector describ- 
ing the contribution of principal components shows a very strong 
contribution from the first 2 components with an exponential decay 
which levels off after the 7th principal component. At this level, 
the detected feature becomes consistent with the shot noise in the 
spectra. The performance of NIR detectors is clearly of great im- 
portance if we wish to detect small signals dominated by a nearby 
bright object. More modern detectors such as the 102 4 2 Aladdin III 
array used in conjunction with NIRSP EC at Keck 2 dMcLean et al.l 
1998) do not show the characteristics dMcLean et alj|2000l) of the 
older Aladdin II array used to secure the data in this work. Fu- 
ture de tectors promise even greater stability with Bezawada & Ivesl 
(2006) finding 0.05 per cent flatfield stability over a period of one 
day for a 2048 2 HAWAII- 1RG engineering array. 
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Figure 6. As for Fig.[5]but using simulated data (50 spectra taken at two epochs) with a mean S/N ratio of 300. The dashed sinusoidal curve represents the 
known K p of the simulated planet. The spectra were simulated to mimic the multi-order capability of existing cross-dispersed spectrographs, covering 7 orders 
and spanning the region of 1 .90 |am to 2.45 |am. Resolutions of R ~20000 (top) and R ~40000 (bottom) were used. A detection at a level of ~ 3 times 
the 99.9 per cent confidence level is achieved, indicating a limiting detection (with the same confidence) for spectra obtained with a mean S/N ~ 170. The 
greyscale in the \ 2 plots runs from 0.99 (white) to 1.00 (black). 



Current observational evidence for a planetary signal at 
2.2 am therefore remains inconclusive. Initial attempts to search 
for extrasolar planetary spectra (by the method of difference anal- 
ysis where the spectrum of the host star seen during eclipse of 
the planet is subtracted from the spectrum of the combined star 



an d planet, seen out of ecli pse) revealed null results. The analysis 
by iRichardson et alj |2003) failed to detect the predicted contin- 
uum 'bump' a t 2.1 -2.2 p.m at a level of logio (F p /F t ) — —3.52. 
ISnellerJ Eowb did not detect the eclipse shape from 2.3 (xm sec- 
ondary eclipse observations of HD 209458b, but nevertheless the 
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mean drop in flux (based on two measurements) of 0.115 ± 0.139 
per cent indicates a planet-star flux ratio of 0.0010 ± 0.0009. 
While this confirms model predictions for a planet-star flux ratio, 
iDeming et al.l d2006h note that the size of the uncertainty does not 
exclude the blackbody result. 

While advances have been made in detec ting spectral fea- 
tures in the mid-infrared 7 .5 - 13.2|j.m interval JRichardson et al.l 
l2007l : lGriHmair et alj2007h using Spitzer space telescope observa- 
tions, models would suggest that the higher contrast ratio in the 
NIR will necessitate the kind of signal enhancement techniques 
used in this work. It is clear that ther e is still much t o learn 
abo ut the atmospheres of CEGPs. Both [Richardson et al.l d2007h 
and lGrillmair et all d2007l) have found that there is little evidence 
for the H2O opacities present in model spectra. Ric hardson et al.l 
d2007h however find evidence for silicate clouds in the spectrum of 
HD 209458b. 

The simulations we have carried out in this work make the 
assumption that the opacities in current models are correct, even if 
their strengths are not. If many species are omitted or are simply not 
present in reality, this could have significant impa ct on the detection 
metho d, either for better or worse. As with the Richardson et al. 
(2003) search for the 2.1 -2.2 11m bump, our models contain unre- 
liable OH opacities and oscillator stren gths and missing CH4 opaci- 
ties which are seen in T dwarf spectra (Bur gasser et al. 2006). This 
is clearly important when using methods which aim to detect the 
'continuum' shape, while we found (Sj5]l that incorrect positioning 
of greater than 15 per cent of lines would affect our ability to make 
a clear detection. Conversely, competing models could be used to 
obtain an optimal extraction of the absorption signature of a plan- 
etary spectrum. Switching off model opacities for species which 
are known to have unreliable positions may be necessary to en- 
able recovery of the planetary signal. In the case of a clear detec- 
tion, it should then also be possible to distinguish between differing 
phase function models by monitoring the improvement in \ 2 m our 
matched filter analysis. 

If the albedo spectrum is largely independent of wavelength in 
the regions where absorbing spe cies do not play a sig nificant role, 
including 2.2 |i.m (see Fig. 8 of lSudarskv etal. (2000)), we would 
expect a very small reflected light component in the NIR, given 
the upper limits found in the optical. Our data contain only 10 
lines arising from the spectrum of HD 75289 itself, with a mean 
depth 2.7 times the mean depth of the 98 lines we expect to 
see in the planetary atmosphere. Nevertheless, if the 10 lines are 
seen as a component reflected from possible clouds in the atmo- 
sphere of HD 75289b, we expect little gain in the signal if the 
upper albedo limit of 0.12 iLeighetal.ll2003rJ) is taken into ac- 
count. The maximum reflected light signal in this case is only 
logio(F p /F») = —4.4, over an order of magnitude less than the 
NIR logio(i ? p /i ? ») ~ —3.0 estimate adopted throughout this 
work. 

6.1 Summary 

We have presented an analysis and feasibility study for the detec- 
tion of the NIR spectroscopic signature of the close orbiting extra- 
solar giant planet, HD 75289b. We do not detect the planet with the 
present Gemini/Phoenix data but through simulations based upon 
models and observational constraints we expect t hat current cross- 
dispe rsed spectrographs, such as IRCS/Subaru {Kobavashi et al.l 
l200d) and NIRSPEC/Keck 2 dMcLean et alj|l998h . operating at 
resolutions in the region of R ~ 20000 are sufficient to obtain a de- 
tection with 99.9 per cent confidence. Facilities which could offer 



higher resolution and broad wavelength coverage will be an enor- 
mous benefit. With a large spectral capability, it will be possible 
to determine the NIR spectral energy distribution in J, H & K re- 
gions, complementing observations currently being made at longer 
wavelengths. 
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